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Abstract: We prove the existence of a fully nonlinear consei^ved curvature perturbation on large scales in 
Galileon-type scalar- field models in two approaches. The first approach is based on the consei^vation of energy- 
momentum tensor of the Galileon field, which is also the familiar approach in understanding the conservation 
in /c-essence or perfect fluid models. We show that the fluid corresponding to the Galileon field becomes perfect 
and barotropic on large scales, which is responsible to the conservation. The difference from A;-essence model is 
that, besides the energy-momentum conservation, the Einstein equation must be employed in order to complete 
the proof of barotropy. In the second approach, we derive the fully non-perturbative action for the curvature 
perturbation ( in Galileon models on lai^ge scales, and argue that ( = const is indeed an exact solution on lai^ge 
scales. This conservation of curvature perturbation is important since it relates the later and the primordial 
universe. 
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1 Introduction and motivation 

Inflation [1] is one of the most exciting and successful ideas in modern cosmology. Over the years, inflation has 
become so popular because of its prediction of a nearly scale-invariant primordial density perturbation, which 
was generated and frozen during inflation to seed wrinkles in the Cosmic Microwave Background(CMB) [2] 
and today's Large-scale Structure (LSS). 

The most popular- and widely studied inflationary models are mostly based on scalar fields [3]. Among 
them, /c-essencial scalar models [4] have attracted much attention due to their broad generality, since which can 
be viewed as the most general scalar- field theories with Lagrangians containing derivatives up to the first order: 
C = C{4>, V(/>). Indeed, this is the "normal" construction of field theories. Nevertheless, one can make a step 
further to consider higher order derivative theories. However, these theories are often suffered with unexpected 
degree(s) of freedom such as ghost(s), breakdown of causality as well as instability problems. This fact makes 
the construction of scalar field models beyond the first-order derivative be not trivial. 

In the study of DGP theory [9] as well as some consistent theories of massive gravity [10], a temi (V(/>)^n(/> 
arises in the decoupling Umit, which yields a second-order equation of motion and thus prevents possible extra 
degree of freedom and ghost. This is the simplest example of a higher-order derivative scalar field Lagrangian 
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which still keeps the equation of motion as second-order (we refer to this property as "second-order preserv- 
ing"). Inspired by this result, the so-called "Galileon" was proposed (and also named) in [16] in flat background 
and extended in [17] in general background, which introduced more general "second-order preserving" terms 
which are nonlinear in second-order derivatives. 

Historically, the most general scalar-tensor theories with second-order equations of motion in four dimen- 
sions was firstly explored by Hordeski in 1974 [18]. This idea — keeping the equation of motion as second order 
— was further developed independently and extended in [19] following the concept of "Galileon", in which the 
most general scalar field theories which we refer to as the "generalized Galileons", whose Lagrangian contains 
derivatives up to the second order £ = V(^, VV(^), but still keeps the equations of motion as second order 
as well as linear in the second-order temporal derivative (j) in arbitrary dimensions, was constructed. It can be 
shown that Homdeski's theory is included in and equivalent to the "generalized Galileons" [19, 23]. Special 
cases of the "generalized Galileons" can also be realized as probing brane embedded in a higher dimensional 
spacetime [20]. The "generahzed Galileons", opens up new windows in the study of cosmology and inflation- 
ary model construction, which have not been explored before. Since the invention of Galileons, a number of 
applications have been made to various contexts in cosmology [21-23]. 

Before applying Galileon models to various aspects of cosmology, however, there is an important issue 
should be investigated: the existence (or non-existence) of a conserved curvature perturbation on super-Hubble 
scales. As in the mostly accepted picture of cosmology, the evolution of cosmological perturbations can be 
divided into three stages: 1) generation during the early period of inflation as quantum fluctuations, 2) exiting 
the Hubble radius during inflation and becoming super-Hubble classical perturbations, 3) finally re-entering 
the Hubble radius in the later radiation/matter dominated era. CMB and LSS formed in the third stage, in 
which we need to set the "initial conditions" for their evolution. These initial conditions are typically given by 
simply the primordial information, which is evaluated in the first stage. What ensures this procedure is that, the 
curvature perturbation is conserved on super-Hubble scales, which allows us to relate the primordial and the 
later perturbations directly. 

The previous understanding of conservation of curvature perturbation highly relies on the study of perfect 
fluid or fc-essencial scalar field. It is well-known that for perfect fluid, there exists a fully nonlinear curvature 
perturbation, which is conserved on large scales to fully nonlinear perturbative orders, if the pressure is a func- 
tion of the energy density, i.e. if the fluid is barotropic [8, 13]. Equivalently, the corresponding perturbation is 
called adiabatic. On the other hand, for perfect fluid or fc-essencial scalar field, scalar perturbation always be- 
comes adiabatic on lai^ge scales [31-33]. This fact explains the existence of a conserved curvature perturbation 
in such models. 

From the point of view of /c -essence [4], the Galileon field is quite strange: it is a model of a single scalar 
field, whereas its energy-momentum tensor takes the form of an imperfect fluid [25, 26]. This is contrary to 
the case of fc-essence scalar, whose energy-momentum tensor can always be cast into a perfect fluid form in 
its own comoving frame. Furthermore, as we will show, the energy-momentum tensor of Galileon field has 
complicated dependence on both scalar field quantities (f), (p etc, but also on gravitational quantities such as 
Riemannian tensors. This reveals that the energy-momentum tensor is not simply controlled by the property of 
the Galileon field itself, but also depends on the whole dynamics of the gravity-scalar field system due to the 
"kinetic coupling (braiding)" between gravity and scalar field, as was also pointed in [25, 26]. This fact makes 
the equation of state, i.e. the relation between p and P much involved. Combining the above complexities 
together, one may wonder if the conservation law for the curvature perturbation still holds in Galileon models. 

'Actually, General Relativity is another elegant example, since the Ricci scalar R contains second order derivatives of the metric, but 
the resulting Einstein tensor in the equation of motion is also up to the second order derivatives. This is because the second derivatives 
of the metric enter R linearly. 
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The present work is devoted to answer this question. We will show the presence of a conserved non- 
linear^ curvature perturbation on large scales in generalized Galileon scalar field models [19]. Two different 
approaches are taken to prove this. The first approach is based on the energy-momentum conservation, which 
is the familiar approach in understanding the conservation in fc-essence or perfect fluid models [8, 13]. We will 
employ the "covariant approach" (or "geometric approach") to cosmological perturbations [13] (see [11, 12] for 
reviews and [14] for related recent development), which can be viewed as a "middle- way" between the usual 
"coordinate-based perturbative" approach (see [5] for comprehensive reviews) and the "gradient expansion" ap- 
proach [6-8]. Both the "gradient expansion approach" and "covaiiant approach" can produce non-perturbative 
results. As mentioned before, an important conclusion got from the previous analysis is that, the existence of 
a conserved curvature perturbation is the consequence of the consei^vation of energy-momentum tensor [8, 13]. 
Precisely, there is a conserved curvature perturbation on large scales and on all peturbative orders, as long as 
the "non-adiabatic pressure" and the "dissipative pressure" vanish, since which act as source terms in the evo- 
lution equation for the curvature perturbation. Correspondingly, the vanishing of dissipative pressure requires 
the fluid to be perfect, and the vanishing of non-adiabatic pressure implies the fluid is barotropic. As we will 
show, this is just what happens for the fluid corresponding the Galileon field on large scales. To this end, we 
divide our proof into three parts: 

• First we show the "Galileon-fluid" is perfect on large scales, or more precisely, the dissipative pressure 
is vanishing. As we have mentioned, there is always a comoving frame for /c-essencial scalar field 
in which its energy-momentum tensor takes the perfect fluid form. This is not possible for Galileon 
field. However, thanks to the higher-order derivative construction of Galileon model, both the energy- 
flow and the anisotropic pressure are suppressed by spatial derivatives. Consequently, the corresponding 
dissipative pressure is suppressed by spatial derivatives just as in fc-essence model. 

• Then we show the "Galileon-fluid" is controlled by two variables (p and (j) on lai^ge scales, which is 
essential in our whole analysis. This is explicit in the A;-essence model, since which depends on (p and its 
kinetic term X = — (V(/))^/2 and thus c/) on large scales. The situation is much more involved in Galileon 
model. As we have mentioned, its energy-momentum tensor depends not only on 0, (j), but also explicitly 
on as well as various gravitational quantities. However, after employing the Einstein equation as well 
as the energy-momentum consei^vation, we can re-express all the other quantities in terms of (j) and (p. 
Essentially, this is because the Galileon model has only one dynamical degree of freedom, the scalar field 

• Finally, we show the entropy/non-adiabatic perturbation^ in the Galileon model is suppressed on large 
scales, as what happens in A; -essence or perfect fluid models [31-33]. Correspondingly, the non-adiabatic 
pressure vanishes and only the adiabatic mode survives on large scales, which implies the conservation 
of curvature perturbation. 

These three steps complete our proof. This approach — relying on the energy-momentum tensor and thus the 
fluid picture of the Galileon field (see [26] for a recent investigation of the fluid picture of the Galileon term 
G{X, (/))□(/)) — is physically transparent and meaningful but mathematically complicated. Thus we develop the 
other approach, that is to derive the evolution equation for the curvature perturbation ( on large scales directly 
and to see if it possesses a constant solution. This is more familiar to the usual coordinate-based approach. 

In [24], the conservation of curvature perturbation in Galileon model with (/>)□(/) term in the La- 
grangian is studied, where the key assumption (p = (/){(/)) is made in order to prove the conservation. We wish 

^Here we refer to the curvature perturbation as"nonlinear" simply mean it is conserved on all orders in a perturbative analysis, 
''since the Galileon model involves only a single scalar field, this is often referred to as the intrinsic entropy/non-adiabatic pertur- 
bation. 
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to emphasize is that, as we show in Sec. 2 and in Appendix A, it is the suppression of entropy mode on large 
scales, which acts as a source term for the evolution of curvature perturbation, that ensures the conservation of 
curvature perturbation. Of course, as a result, the scalar- field is driven to the attractor with cj) = on lai^ge 
scales. 



2 Adiabaticity on large scales 

2.1 Energy-momentum conservation and evolution of curvature perturbation 

We employ the covariant approach to cosmological perturbations [13], and briefly collect the relevant formulae 
below (see [11, 12] for reviews). The covariant formalism is powerful in that, it provides an explicit and 
especially a covariant separation of an object into its timelike and spacelike parts, which is most convenient for 
our purpose for investigating the large-scale behavior of a system, where we will concentrate on the leading 
terms in the expansion with respect to the spatial derivatives'* . 

The covariant approach to cosmological perturbations is highly related to the so-called "3-1-1" or covariant 
fonnalism [15], which splits all covariant quantities into temporal and spatial parts by choosing a special fo- 
liation structure of spacetime. This splitting is usually realized by specifying an arbitrary unit timelrke vector 
field Ua- The energy-momentum tensor can be decomposed as^ 

Tab = PUaUb + Phab + 2n(„gb) + TTab, (2.1) 

where hab = gab + "^aUb is the spatial projection tensor which is orthogonal to the fluid velocity Ua, p and P 
are the energy density and pressure respectively, the energy flow qa and the anistropic stress TTab are both spatial 
tensors which satisfy qau"- = 0, -KabU^ = 0, -Kab = T^ba and g'^^nab = h^-^iTab = 0. The fluid variables can be 
read from (2.1) as 

p = U^U^Tab, (2.2) 
P = \h''''Tab, (2.3) 
qa = -hiu^n,, {lA) 

TTab = {Khf - hl^''hai}j T,d (2.5) 

The spatially-projected derivative Da acting on an arbitrary tensor is defined by 

DaT-' = hi---hl,---hi--- Va'T-'' (2.6) 

where the operator Da is indeed a real covariant derivative, but associated with hab instead of gab- We also 
introduce the decomposition of derivative of u", 

DbUa = VbUa + Ubtta = ^&hab + (^ab + ^ab, (2-7) 



''The covariant approach, which relies on the "3-f1 decomposition" — contrary to the most popular "coordinate-based approach", 
which specifies a coordinate system in which the background is spatially homogeneous, breaks the general covariance in a "hard" 
manner and yields the problem of gauge artifacts — splits the time and space in a covariant manner without introducing any prefeiTed 
coordinate system. This fact allows the covariant approach to grasp the physical transparency which is often lost in coordinate-based 
calculations, and to yield many interesting understandings on cosmological perturbation, including the conservation of super-Hubble 
curvature perturbation. However, the main obstacle in covariant approach is that, it is not easy to perform a perturbative calculation, 
which is needed for practical purpose. 

^Throughout this paper, symmetrization is normalized, i.e. A(^aBb) = ^{AaBt + AtBa). We use a,b, - ■ ■ as covariant indices, 
which are popular in the literatures on covariant formalism. 
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where Q = V is the expansion, = u^VbUa is the acceleration vector, aab = CTba is the symmetric 
tracefree shear tensor, ooah = —^ha is the antisymmetric vorticity tensor. Note the Frobenius' theorem implies 
^ab = if is hypersurface orthogonal, which is just the case we are interested in. 

The projection of the consei^vation of energy-momentum tensor onto u°', i.e. Ua^bT"'^ = 0, yields a 
continuity equation for the energy density and pressure: 



p + e(p + p) = p, (2.8) 



where the dissipative term reads 



V = -q'^aa - Daq'' - n^'aab- (2.9) 



In this section, an over dot " " denotes the Lie derivative with respect to u"". £u, which has the natural 
explanation as a "covariant" temporal derivative. When acting on scalar quantities, it is just / = £uf = 
u"-Vaf- We can define the covector^ [11, 13] 

a 1 
Ca = Daa--DaP, with d = -G, (2.10) 

P 3 

which can be explained as the curvature perturbation. It was shown in [1 1, 13] that (2.8) can induce an evolution 
equation for (a, 

Ca = + S,) , (2.11) 

3p 

where 

r„ = DaP - -DaP, (2.12) 

p 

and 

Sa = (^DaV - jDaP^ + ^ (^DaO - -DaP ] , (2.13) 

are the nonlinear non-adiabatic pressure and dissipative pressure respectively^ . Now it becomes explicit from 
(2.11) that, in order to prove the conservation of (a, i-C- to show it satisfies a conservation law (a = 0, the 
source term in the right-hand-side of (2.1 1) must vanish, i.e.: 

ra + S„ = 0. (2.14) 

As we will see, this is just what happens in Galileon model on large-scales, where both the dissipative pressure 
Sa and the non-adiabatic pressure Fa vanish sepai^ately, which implies the fluid corresponding to the Galileon 
field becomes perfect and barotropic on large scales. Correspondingly, the scalar perturbation in Galileon 
model becomes adiabatic. 




2.2 Galileon model with G{X, < 

In this section, we make a detailed investigation on the Galileon model with G{X, This term has already 

included the essential ingredients in our analysis, but is simple enough in order to help us to grasp the physical 



*In the covariant approach, perturbations are usually defined as spatial gradients, e.g. the scalar perturbation Scf) in coordinate-based 
approach is represented by Da(l>, the energy density perturbation Sp is represented DaP, etc. See [11, 13, 15] for details. 

^Throughout this paper, we use the words "dissipative pressure" or "entropy perturbation" to indicate the effects of imperfect fluid 
nature of the energy-momentum tensor of the "generalized Galileons", which is usually used in the cosmological context. One should 
keep in mind that, the Galileon model is a closed and essentially conservative system, which implies there is no real physical dissipation 
or entropy production. 
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picture without being immersed in cumbersome mathematics. A full treatment of the generahzed Gahleon is 
presented in Sec. 2.3. 

Consider the Galileon model 

C^ = K{X,cl)) + G{X,(l))U(l), (2.15) 

where X = — (Va(/>)^/2. The con^esponding energy-momentum tensor Tab = ~"^7=~^^^^^^ of the scalar 
field reads 

Tab = Agab + BVa(j)Vb(j) + VaCpVbG + VbCpVaG, (2. 16) 

with 

A = K -Vafpy^'G, (2.17) 
B = K,x + G,xa(j). (2.18) 

Here and in what follows, K^x denotes dK/dX, etc. (2.16) does not takes the form of the energy-momentum 
tensor of a perfect fluid, due to the presence of G{X, (p). In the covariant formalism, the corresponding fluid 
quantities can be calculated sti^aightforwardly: 



A + B(P^ + 2G,x(l)X, (2.19) 

(2.20) 

B(l)DaCb-Gy{XDa<P + <PDaX], (2.21) 



P = A + ^ [BDa(t)D''4> + 2G,xDa(pD''X 



TTab = BDaCpDb<t) + G,X {Da<l)DbX + Db^DaX) 

[BD,4>D'(t) + 2G^xD,4>D'x) hab, (2.22) 

where B = B + 2G^^, f = u'^V af- (2.19)-(2.22) can match the corresponding results derived in [26] in 
coordinate-based approach. 

It is important to note, as in the single /c-essencial scalar field model, both the energy-flow Qa and anisotropic 
pressure -Kab involve at least one spatial derivative, and thus can be neglected on large scales. Precisely, we may 
write 

qa ~ 0{D), TTab ~ 0{D^), (2.23) 

where 0{D^) denotes n-th order in spatial derivatives. One difference between Galileon and fc-essence is that, 
in single field A;-essence model, we are always able to choose an unique "comoving frame" for the scalar field 
defined by Da(l) = 0, in which the fluid becomes exactly a perfect one, i.e. = tt^™™^ = 0. This is not the 

case for Galileon. As from the expression for qa and TTab, in the comoving frame for the scalar field, tt^™""^ = 
as in /c-essence model, however the energy-flow^ 

^(com) ^ _G,x4>DaX / 0, (2.24) 

due to the dependence of G on X. As was firstly pointed out in [26] and will be shown in the following 
sections, this is also the case in the full "generahzed Galileons". This fact prevents the energy-momentum 
tensor of Galileon from a perfect one. 



*This phenomenon was also found in [22] in studying the scalar perturbations of Galileon model (2.15), where it was firstly pointed 
out that, contrary to the fc-essence model, the "comoving gauge" Tf = and the "uniform field gauge" Stj) — are not equivalent in 
Galileon model (2.15). 
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For our purpose in this work, we concentrate on the energy density p and the pressureP, whose explicit 
expressions are 



p = p{4>,4>,Q,A^p^) = -K + i^K^x + G,^ - G,x<pe) + V^p), (2.25) 
P = P{<j), 0, 4>, ^(P)) = K+ (g,^ + G,x4) + ^(P), (2.26) 



where 'D(p) and 'C(p) are terms at least second-order in spatial derivatives, i.e. I?(p) ~ 0{D'^) and ^^(p) ~ 
0{D^), which can be safely neglected in the large-scale approximation, whose expUcit expressions can be 
found in (C.1)-(C.2). In (2.25)-(2.26), the arguments Af^p-^ and ^(p) formally remind us the dependence of p 
and P on all other terms which are higher-order in spatial derivatives. As was pointed in [25, 26], besides 
those terms which are higher-order in spatial derivatives, the energy density and pressure have the following 
functional form 

p = p{<p,4>,e), p = p{<p,4>J)- (2.27) 

It is useful to recall that in A;-essence model, both the energy density and the pressure take the form p = p{X, (/)) 
and P = P{X, (f), which can also be written as p = p{4>, (p) and P = P{4>, cj)) up to terms which are higher- 
order in spatial derivatives. In Galileon model (2.15), due to the dependence of P on cj), the pressure P is 
determined only after solving (p in terms of cj), cj) etc, by making use of the dynamical equations of motion. The 
dependence of p on Q, however, is essential for the Galileon model, which implies the energy density depends 
on the dynamics of the full system, rather than solely on the property of the scalar field. For later convenience, 
note in(2.25)-(2.26), p is linear in B and P is linear in (p with coefficients: 

P,e = ;:^7;; — -G^x(p , (2.2«) 



de 



dP(<p,(f>,(l>,A(p)) .„ 
Pi ^ I : ^ ^ = G^xcP^ (2.29) 



both of which are proportional to G^x ■ 

From (2.23), Qa C>{D) and i^ab ~ 0{D^), the definition of the dissipative term (2.9) implies V ~ 
0{D'^), and thus in Galileon model (2.15), the dissipative pressure is third order in spatial derivative^, 

(2.30) 

and thus can be safely neglected on large scales, which is similar to the case in (even multi-field) /c-essence 
model [13]. Thus in order to show the conservation of Ca on lai^ge scales, it is essential to show Ta ~ 
on lai^ge scales, which implies the fluid corresponding to Galileon (2.15) becomes barotropic on large scales, 
equivalently, the scalar perturbation becomes adiabatic on large scales. 

In the appendix A, we discuss the condition for the barotropy of the "2-2"-mapping system with p = 
/?((/>, 0) and P = P{(p,(l)), which are explicit for /c-essence model. Here we will show that, if the universe 
is dominated by Galileon field, the energy density and pressure of Galileon can always be cast into the "2-2 
mapping" form, which is essential in our analysis. To this end, we have to make use of the Einstein equation 
Gab = Tab- Precisely, by projecting both sides on u°'u^ we arrive at the energy constraint 

u'^u'Gab = p. (2.31) 



'Readers who are familiar with the "3+1 formalism" would recall Ua = Da InN, where A'^ is the scalar lapse function. Thus the 
acceleration aa is first order in spatial derivative. 
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Using the decomposition (in writing this, we assume u°- is hypersurface orthogonal, which is actually the case 
of cosmological interest) 

DhUa = -Qhab + (Tab, (2.32) 

as well as the Gauss-Codacci equations (see Appendix F), we are able to rewrite the energy constraint as 



with 



n=l{-'R- UahO 



ah 



(2.33) 



where is the intrinsic Ricci scalar of the spacelike hypersurface orthogonal to the 4-velocity Ua, which is 
just the Ricci scalar associated with hah- (2.33) is the covariant (and non-perturbative) version of the familial" 
background equation ?>H^ = p. In an almost FRW universe, which is the interest of cosmology, R^^^ ~ 0{D^). 
Similai^ly the shear aah also rapidly decreases on large scales. Using (2.33) and (2.25), we can solve as 



e = e( 



.A 



(0). 



n, 



(2.34) 



where "it" denotes different possible choices of branch, ^(e) denotes the dependence of G on all terms which 
are higher-order in spatial derivatives. In (2.34), it is explicit that on large scales, B depends only on (p, 0. 
Using the definition of p together with the energy constraint (2.33), we can eliminate Da@ and 6 and get 

-1 



DaP 



P 



1 



3jO,e 
26 

26 



p,^4> + p,4>^ + pAip)\p) - 



(2.35) 
(2.36) 



Note even on large scales, DaP 7^ p ^Da(j) + p^^Dac/) etc, due to p^e 7^ 0. As we have mentioned before, 
we can solve from the dynamical equation of motion for the salar field, or more conveniently through the 
energy-momentum conservation (2.8) together with the energy constraint (2.33) to get 

902 



= ^(0,0,A(^)) 
with 



(1 + Ar 



PA, 
P.s' 



+ 



1 



A 



@P,e 



3p,e 
26 

26 



A> + 2G,^ - G,x06 + 



(2.37) 



(2.38) 



Pj^^^ is higher-order in spatial gradient which is given in (C.3). Again A^^-^ denotes the dependence of on 
terms with higher order spatial derivatives bisides and (j). Thus we have 



DaP 
P 



Pi + PiA 



Da(l)+ P,, + P2A,s] D, 



+ P:,A,A„DaA^^;j^^ + P^A^^^DaA^p), (2.39) 



PlA Ar- A, 



-0 '^(<*) (0) 



(2.40) 



Although (2.35)-(2.40) look cumbersome, the meaning of whose structures is simple: besides all the terms 
which are higher order in spatial gradients, e.g. ^(p), etc, on large scales, both p and P are only determined 
by and (j), which is essential for our following analysis. 

Having shown that both the energy density and the pressure can be written as functions of and together 
terms which can be neglected on large scales, we are at the point to show that ~ on large scales. For later 
convenience, we introduce the "relative entropy perturbation" between any two variables X and Y: 

DaX DaY 



Sa{X,Y) 



X 



Y 



(2.41) 
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which will significantly simplifies our following discussion. In terms of Sa{X, Y), the non-adiabatic pressure 
is just 

Ta = PSaiP,p). (2.42) 
As usual, we introduce the so-called "comoving density perturbation" 

€a = DaP - iDa<P = pSa{p, 4>). (2.43) 

It is well-known that in the perturbation theory of perfect fluid or A: -essence scalar field model, the comoving 
density perturbation is suppressed on large scales. Actually as we will show, what is suppressed is the following 
quantity [13] 

ea = DaP - Oqa. (2.44) 
To this end, we have to make use of the Einstein equation. The momentum constraint is given by 

u'^GbcK = -Qa, (2.45) 
which can be recast into a more convenient form 

D''aba - ^DaQ = -qa, (2.46) 

through the Gauss-Codazzi relations (see Appendix F). Combining the energy constraint (2.33) and the mo- 
mentum constraint (2.46) together, we get 

DaU + QD'^CTha = DaP - QQa = (2.47) 

which is the covariant generalization of the Poisson equation in linear theory. (2.47) is one of the essential 
equations in our analysis. In an almost FRW background, the left-hand-side of (2.47) is suppressed by spatial 
derivatives, which implies that « on large scales. 

In perfect fluid or fc-essence model, it was shown that ea and ia exactly coincide on large scales (or in 
linear theory) [13]. In Galileon model, as we will show, ea and ia do not coincide but are proportional to each 
other. Explicit manipulations show that (see Appendix D for detailed derivation) 

ia = {l + X)ea + {l + A) + , (2.48) 

where A is given in (2.38), and 'C(g^) are higher-order spatial derivative terms, whose expression are given 
in (D.6) and (D.8). On large scales where ~ ^(?a) ~ we have 

ea-{l + X)ea. (2.49) 

For fc-essence model, p does not depend on the expansion G explicitly, thus p e = and A = 0. Thus (2.49) 
is the generahzation of the relation ia ~ ea in A;-essence model to the case of Galileon (2.15). Similarly the 
non-adiabatic pressure can be evaluated as 



Ta = rCa r2?(,^) + U - ^^7^ -Sa, (2.50) 



where 

C = P,4>P.^ - P.^P<t> + P.6 ( pM.^ - P.6A'I> ) , (2-51) 
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and Sa is higher-order in spatial derivatives and is given in (D.4). On large scales, (2.50) impUes 



(2.52) 



(2.49) and (2.52) explicitly show that on large scales Ta, and ia are proportional to each other. Whereas we 
have proved that on large scales is suppressed (see (2.47)), thus we get the conclusion that, on large scales, 
the non-adiabatic pressure is vanishing 

Ta ~ 0. (2.53) 

This fact ensures the conservation of curvature perturbation on large scales. 

Let us try to understand the above result. The Galileon model (2.15) involves only one dynamical degree 
of freedom, the scalar field (p. Thus on large scales, the system involves only two perturbation modes Dac/) 
and Da4>, which can be combined linearly to give one adiabatic mode and one entropy mode. The fact that the 
system contains only one entropy perturbation implies that. Fa, ia etc, although defined differently, must be 
equivalent characterizations of the same entropy perturbation mode, and must be proportional to each other on 
large scales. What is essential is that, the entropy mode has no dynamics for itself and is suppressed on large 
scales, only the adiabatic mode survives, as what happens in the perfect fluid or fe-essence models. To conclude, 
as in single-field A;-essence model, the scalar perturbation in Galileon model (2.15) becomes adiabatic on large 
scales and the corresponding curvature perturbation (a is conserved. 

2.3 Generalized Galileon model 

The above analysis can be applied to the generalized Galileon model proposed very recently in [19], although 
the calculation becomes dramatically cumbersome due to the increasing complexity in the structure of the 
energy-momentum tensor when including higher order Galileon terms. 
The full Lagrangian of the generalized Galileon is given by [19] 



with 



£(0) ^^(0) 

= i^w (x,0)n,/., 

£(2) = i^(2) (X, 0) [(0^)2 _ (V,Vfc0)2] + R (X, 4>) , 

£(3) ^ ^(3) ^) Jp^)3 _ (V,Vb<^)2 + 2 (VaVbCl))^ 



where X = -(V(/))V2 and 



Q(2,3)(^^^) = dY K^^'^^ {Y, 

Jo 



(2.54) 

(2.55) 
(2.56) 
(2.57) 

(2.58) 

(2.59) 



that is = K^^'' K^'^Hx,(j)) etc are arbitrary functions of X and (p, R and Gab are Ricci scalar and 
Einstein tensor respectively. 

The energy-momentum tensor of the generalized Galileon model (2.54) is given by 



n=0 



(n) 
ab ' 



p(n) 

" ab 



2 



ab 



(2.60) 
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where a superscript "(")" denotes the contribution from We have 

T^^ = K^'^^gab + K%^Va^Vt^, (2.61) 

+ (K?n</' + 2K^l^^ Va(l>Vb(l> + 2K%^V^,(l)Vb)X, (2.62) 

where r^°^ + T^^^^ is just (2.16) with iT^o) = K and = G. Contributions from £(2) and £(3) are much 
involved, which are given by 



-AK^l^VcCpV'X - 200 (k^JVc^V^X - 2XK'^f + qJ^^) + 2K(2)i^^^v^(/.VV} 



+ 



(2) , OZ<'(2)^ Y7, ^T7., ^ , OI^(2) 



Va(/'V6(/. 



+4 (a^K';^> + 2Kl;>) V(,</)V,)X + 2Kl^' {VaXVbX - 2V,XV=V(„ 
+2 (i^(2)Q^ ^ i^(J)Ve<AV^X - 2Xi^^(^') + Q^)) VaVb0 

-2i^(2) f^c^^^^^^^^ ^ 2V(„(/.ii,),V^(/. + iJeadfcV^c/^V^c/.) - 2Q(2)^^,, (2.63) 



(3) (I) (2) (3) 

and T^f, is given in Appendix E. As T^^' depends on 0, T^^' and T^;, have exphcit dependence on gravitational 
objects such as R, Rat and Rabcd etc. 

Now we concentrate on the large scale behavior of the energy- momentum tensor. To this end, one would 
find the following two expressions be useful 

Va<l) = -Ua<j) + DaCp, (2.64) 

VaVb<P = UaUb (0 - a^I^c^) - 2n(^ (d^^^ + D'uh)Dc4)^ - DbUa<j) + DaDb<l), (2.65) 

where DbUa is given by (2.32). By keeping only the leading terms in spatial derivatives, we can write 

Va(t> = -Ua^ + 0{D), (2.66) 
VaVfoc/) = UaUb4> - ^e4>Kb + 0{D), (2.67) 

where again, 0{D) denotes higher-order spatial derivative terms which will be neglected below. For later 
convenience, note 

X = ^^^ + 0{D^), VaX = -Ua<P4> + OiD), (2.68) 

which can be verified easily. At this point, it is useful to note, besides u"" and hab, ^a4' and X depend on (j), 
while V aSJb4' arid depend on (j), (p and on large scales. 
Plugging (2.66)-(2.67) into (2.61)-(2.62), we have 

- (i^SV^ - K('^) naUb + K^'h^b, (2.69) 

- (^S^ - 0<^) ^^^aUb + [k^x^4> + K^l'^) ^^Kb, (2.70) 

from which the corresponding large-scale expressions for the energy density and pressure p^^\ p^^\ P^^^ and 
p(i) 

can be read easily. From (2.69)-(2.70), it is explicit that t|||^ and t|^^ become that of a perfect fluid 
separately (and automatically) in the large-scale approximation. Actually, (2.66) and thus (2.69) are exact in 
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the comoving frame of where Dacj) = 0, which yields the well-known conclusion that /c-essence scalar field 
in its own comoving frame corresponds to a perfect fluid. 

(2) (3) 

Similar procedures can be applied to T^^ and . In the following, we perform a detailed analysis on 
t'^ . The large-scale expression for is 

Tab - ^aUb 



+hab 



+2,^2g(2)^ + 2ir3) j'^ _ 2^ A ^2 ^ 2 U + ^e<^) QI') 



2E:(2) (2u^^Rb)X<j)^ + i^eadb^^'n'^'A') - 2Q(2)i?^^. (2.71) 



where we have used the decomposition of R, Rdu'^W^ on large scales, see Appendix F for details, have 
almost taken the fomi of a perfect fluid, except the last line. Instead of trying to show the last line in (2.71) 

(2) 

can also be written in a perfect form, it is easier, however, to evaluate the con^esponding energy-flow qa and 

f2) 

the anisotropic stress tt^^ directly and to see if they are indeed vanishing on large scales. Only the last Une in 

(2) (2") 

(2.71) contributes to qa and tt^^ , thus we have 



-Q^^^^hlu'^Rbc + OiD), {1.11) 
-2 {hlhi - hi^''hai}j [K^^^^Recfdu'u^ + Q^^^Rcd) + 0(0^). (2.73) 



Contrary to T^^^ and Tj^l\ at this point it is not explicit if qa ~ TTab ~ and thus t[^^ takes the form of a 
perfect fluid under the large-scale approximation. Now we need to employ the Gauss-Codazzi relations. (F.2) 
yields 

h''XRbc = DbC7',-'^Dae^O{D), (2.74) 

which implies 

gi^) ^ OiD), (2.75) 

and thus vanishes in the large-scale approximation. It is interesting to note, a special case is when the coefficient 
in (2.72) vanishes, i.e. 

i^(2)^2_g{2)^0^ 

In this case q^^ vanishes automatically without using the Codazzi relation. This is satisfied when K^'^^ oc X, 

(2) 

which corresponds the original Galileon model studied in [16]. Similarly analysis can be applied to vr^^ . To 
this end, we need to use the Ricci equation (F.3) as well as (F.4), (F.5) and (F.6), which yield 

hlht - ^h'^'hab^ RecfdU'u^ = (e^ + 36) hab + hab + ©) + ^(^')' (2-76) 



hlht - -h'^habj Red = 3 (0 + ®) Kb - 3 (^36^ + 26 - — - 6J + 0{D'){1.11) 

which cancel exactly and leave us terms of order 0{D'^). Thus, we can conclude that 

TriJ ~ 0{D^), (US) 

as what happens for vr^j^^. Similar- procedures can be applied to xj^f^ and yield the same conclusion. To conclude, 
the energy-flow qa and dissipative pressure TTab in the generalized Galileon model (2.54) ai'e higher-order terms 
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in spatial derivatives and thus can be safely neglected in the large-scale approximation, which implies the fluid 
corresponding to the Galileon model (2.54) becomes a perfect one on large scales. One may notice that up to 
now, the Einstein equation has never been used. Thus the above conclusion is generic, which is irrelevant to the 
underlying theory of gravitation. 

Having shown both T^^ and T^^ become the form of a perfect fluid on large scales, which implies the 
corresponding contributions to the dissipative pressure vanish w 0, the next step is to show the fluids 
described by T^^ and T^^ we, actually barotropic, and thus the con^esponding non-adiabatic pressure also 
vanishes 0. From (2.71), the energy density and pressure from can be evaluated straightforwardly^" 

« -|g(2)G2 + \k^^)q^^ + Ik^'^Q^' - 2Q;><^ - 2J^g)G<^^ (2.79) 



p(2) ~ Iqm (q-^ + 20] - -K^^^ (Q^<i) + 2Q<j) + 26^) - -K%^Q<j)' 



and 



+2Q;S'^' + 2Q?^^ + l^'P^ + 2i^^^'V' (^4> -1@<P], (2.80) 

(3) 

where again we used various relations in Appendix F. Similarly, for T^j, , 

p(3) ^ _}^k^3)q3^3 _ lK%^e^' + GQ^^^Q^'' + 2Kfe^\ (2.81) 



and 

9 V / 3 ' 

-^Q^i©'/'' - 2g;jV (©'0 + 2G</. + 4e(^) + ^Kfe<}>^ (q<P - 64>) ■ (2.82) 

Now we concentrate on the functional structure of p and P in generalized Galileon model (2.54). On large 
scales, note JT^"), etc are functions of ^ and X ^ <j>'^/2 and thus (see (2.69)-(2.70), (2.79)-(2.82)) the 
energy density p and the pressure P take the following general functional form 

p = pi + ep2 + e^P3 + o{D), (2.83) 

P = Pi + ePa + ©^-Ps + 4>iPi + ©^5 + ©^i^e) + 0(i^7 + QPs) + O(^), (2.84) 

where pi and P/ are functions of (/> and whose explicit expressions can be read from (2.69)-(2.70), (2.79)- 
(2.82). Let us focus on the dependence on 0, Q and (j) of p and P. On one hand, p is simply a binomial of 0, 
which does not depend on Q and (p. On the other hand, and cf) enter P only linearly. The point is, this type 
of functional fomis of p and P allows us to solve 0, and (p in terms of (p and (p. To this end, the Einstein 
equation and the energy-momentum constraint (which is equivalent to the evolution equation for the Galileon 
field) have to be employed. Precisely, we get two equations from the Einstein equation: 

p = u'^u'C,, = l®' + l (^^'^ - ^«''^''') = I®' + ^(^)' (2.85) 
3P = h''''Gab = -©^ - 20 - ^P(3) _ laabd"^ + 2aaa'' + 2Daa°- w -0^ - 20 + ©(L'^), (2.86) 



"'Actually, these large-scale expressions for p and P can also be derived through simply the background equations of motion by 
identifying O = 3-ff. 
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where we used relations in Appendix F. Alternatively, one may employ the energy-momentum constraint (2.8). 
In any case, we can solve and (j) in terms of (f) and (j) and thus finally arrive at: 

p = p{4>,4>), P = P{(P,<i)), (2.87) 

just as what we did for G{X, (/))n0 model. It is this fact that both the energy density and pressure can be cast 
into a "2-2" mapping that ensures the barotropy on large scales, as is discussed in Appendix A. The key point 
is that, as what happens for K{X, (p) + G{X, (f))\I\(f) model (see (2.47)), the non-adiabatic/entropy perturbation 
mode is always suppressed on large scales, leaving us only the adiabatic mode. Actually as emphasized in 
Appendix A and B, this is a general conclusion as long as p and P are determined by two variables (p and (j), 
which is irrelevant to the details of the model. 

To end this section, we would Uke to make a general comment in order to understand why we can always 
get p = p{(/),(/)) and P = P(0, 0), which is essential in our analysis. The Galileon action (2.54) and the 
corresponding energy-momentum tensor Tab are generally constituted by the following ingredients: 

0, Va(/>, VaVb4>, Rabcd, 

as well as the metric gab- In the scalar field sector, from (2.66) and (2.67), on large scales Va(f> is proportional 
to (p and V a^h<P is function of cp, cp and Q. In the gravity sector, (F. 1)-(F.6), in all contractions of Riemann 
tensor with VaCp, ^a'^b'P and gab (and thus with hab and u^) which will contribute to the energy-momentum 
tensor, only terms proportional to 0^ and survive on large scales. Thus, the system is apparently controlled 

by 

(p, <p, <p, 0, 0, 

on large scales. However, the Galileon model (2.54) is constructed in such a very special manner that the 
resulting equations of motion are second order, both for the Galileon field and the metric. Two implications are 
in order: 

• Mathematically, repeated coordinate indices only appear once, which implies that the second order tem- 
poral derivative (p enters the theory only linearly^^ , which can be verified explicitly, as 

U(P^ -4)- e<p, (2.88) 
{Dcpf - {VaVb^f ~ 2@<P4> + (2.89) 

[Dcpf - SDcP {VaVb^f + 2 {VaVb^f « -2e^^ - ^Q^cj)^, (2.90) 

all of which are linear in (p. This is also what happens in the energy-momentum tensor as well as the 
equations of motion. 

• On the other hand, the second-order-preserving construction prevents Galileon from any extra dynamical 
degree of freedom, which means is not a dynamical degree of freedom and can always be solved in 
terms of the real dynamical degrees of freedom. Mathematically, this fact implies in the equations of 
motion (and thus in the energy-momentum tensor and Einstein tensor), there is no temporal derivatives 
of higher than the first order; and only appear linearly, as we have seen. 

These two facts allow us to solve both (and thus 0) and cp> in terms of (p and <p, at least in principle. Physically, 
this is simply because the real dynamical degrees of freedom in Galileon model (2.54) on large scales are (p and 

4>. 

"it is also this fact makes the Galileon model has a well-defined Cauchy problem. 
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To conclude, the fluid corresponding to the Galileon model (2.54) becomes both perfect and barotropic on 
large scales. This fact ensures the adiabaticity of the scalar perturbations in our model on large scales and thus 
the consei^vation of curvature perturbation. 

3 Shift-symmetry on large scales 

In the above section, we have shown the existence of a conserved quantity which has the physical meaning as 
the curvature perturbation, using the conservation of energy-momentum tensor of the scalar field. The proof is 
physically transparent and meaningful but is mathematically involved. Thus it would be interesting to see the 
behavior of the curvature perturbation on large scales by deriving its own evolution equation, and to verify if it 
indeed possesses a constant solution. 

One may recall that in A;-essence models, the quadratic order action for the curvature perturbation ( (de- 
fined by 6gij = a^(l + 2()5ij) is 5(2) = / drjd^xa'^^ [('^ — c^{di()'^]. On large-scales (neglecting spatial 
derivatives) the action reduces to an action for one-dimensional system S(2) oc f drja^-^Ch ^'^^ Cl = Cl(?/)> 
where Cl(^) is the long wavelength component of the curvature perturbation, which can be viewed as the av- 
eraged curvature perturbation of a given local Hubble volume. The corresponding equation of motion for is 
(cL^-^cCj ~ 0> which has a dominant solution Cl = const, during inflation. This indicates the conservation 
of Cl, i-S- the conservation of C on large-scales at linear order At this point, it is interesting to note that the 
absence of "mass term" ('^ in its action ensures the conservation of (. Actually, due to the absence of "mass 
term", the large-scale Lagrangian a^-^C^ has the "shift symmetry" C — ^ C + const, which is responsible to the 
conservation of C- 

Following this logic, now our task is to show (^l = const, is actually a solution (dominant solution in 
an expanding background) of the Galileon models on non-linear orders. As we will see, the shift symmetry 
C ^ C + const is preserved on fully non-linear orders, which is essential for the conservation of C- 

We work in uniform field gauge with 6(j) = 0, and since we will study the perturbation on the large scales, 
we take the metric perturbation as 

ds^ = a? (^-e^'^drj^ + e^^dijdx'dx^^ , (3.1) 

where a = a{rj) is the scale factor and C will be identified as the curvature perturbation. The vector modes have 
no dynamics for itself (do not propagate), which decay on large scales in scalar field models. Tensor modes will 
couple to scalar mode on nonlinear orders. However, these couplings arise only through the spatial derivatives 
of scalar modes, e.g. diC, which can be safely neglected for our purpose. In the following, we neglect spatial 
dependence of a and which means we focus on their long wavelength components q;l and Cl> which are 
spatially uniform in a local Hubble volume. In the following we neglect the subscript "L" for short. 

As a warm-up excise, let us consider the simplest case with = X. The background equations are 
QT-L^ = (p'^ and Ti' + = 0, where a prime denotes derivative with respect to comoving time ry. In this case, 
the large-scale Lagrangian is 

C « -3a2e-"+2'^ (C'^ + 2nC) ■ (3.2) 
Vaiying (3.2) with respect to a yields an equation for 

e^^C' i^n + C) =0, (3.3) 

which implies C' = ^i^d thus the conservation of (. At this point, we wish to emphasize that (3.3) is a fully 
non-perturbative result, and thus ( = const, is a non-perturbative solution of the system on large-scales. For 

'^In this section, Lagrangian £, includes \J—g. 
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completeness, it is interesting to solve the constraint a in terms of C- Varying (3.2) with respect to C yields an 
equation: 

d (c' + ^) - md - C - ^c'^ = 0, (3.4) 

from which we can solve 

"='"('4)4/*'cT«- 

Expanding (3.5) perturbatively gives the familiar results [30]: a = ^ + • • • . 
3.1 Galileon model with G{X, (t>)a(f> 

Now we turn to the Gahleon model (2.15). After some manipulations, we arrive at a fully non-perturbative 
action for the curvature perturbation Q and a on the large scales: 



S = j dr]a^e-''+^^ [-3{'H + C'f + a^e^^P (X, 0) + G' (X, , 



(3.6) 



where on large scales X ^ We emphasize that (3.6) is exact in the limit of large scales, by neglecting 

spatial gradients. No other approximation is made in deriving (3.6). Mathematically, as we have already 
emphasized, the large scale action (3.6) only involves C, through an overall exponential factor. Thus under the 
constant shift of C,: 

C ^ C + c, (3.7) 
where c is an arbitrary non-vanishing constant, the Lagrangian changes as 

C e^^C, (3.8) 

which has no dynamical meaning and implies the system is actually invariant under the constant shift of C,. This 
is just the case of A;-essence model we mentioned before, while the presence of Galileon-like term G{X, (l)){^4> 
does not change this property. To conclude, if C = is a solution of the system, C, = const ^ must also be a 
solution. 

Although the above argument has akeady completed our proof, one may convince him/herself by deriving 
explicitly the evolution equations. Varying (3.6) with respect to a and C, yields 

Q = 2,{n + C'f + [n + C') + a'e^^P - {Px + G,^) 4l\ (3.9) 

and 

= + C') {n + 3C' - 2a') + a^e^^P + 2 (?^' + C") 



respectively. On the background level, by definition, a = C = is a set of solutions of equations (3.9)-(3.10). 
Due to the absence of C term in the equations, one can exphcitly verify that actually 

a = 0, C = c, (3.11) 

where c is an arbitrary non- vanishing constant, is also a set of exact solutions of the nonlinear equations (3.9)- 
(3.10). This completes our proof that, on large scales, Q is conserved non-perturbatively. Note the equation for 
C (3.10) is linear in (^", which impUes the system has a well-defined Cauchy initial- value problem. 
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3.2 Generalized Galileon model 



The above analysis can be easily extended to the generalized Galileon model (2.54). After some manipulation, 
the exact large-scale Lagrangian takes the following form 



C « e=^^ [A{C', a) + a'B{C', a) + C"C{C', a)] , 



(3.12) 



where A, B, C are functions of C,' and a as well as background quantities such as a and T-L, whose explicit 
expressions are given by 

A = ^Te"^" 



a^e^" (3(H + C')(^ + 2C') + 3?^' + 6Q(2) ((^ + (,'){% + 2C') + ?^') - K^^^ {{2U + 3C' 
-Qa'^e'^'^in + CO (-i^^^V (CV + 0") + 3g(=^) ((3C'(^ + C') + 2?^') 0' + + Q')(t>" 



B = — e 



-5q 



C 



-3a 



(3.13) 

(3.14) 
(3.15) 



We emphasize that in (3.13)-(3.15), functions K^^^ etc are their corresponding large-scale expressions by ne- 
glecting all spatial derivatives, e.g. K^^^ = K^^\^^-^(f)''^ ,4>) etc. Again, from (3.12) it is explicit that C, 
enters the Lagrangian (3.12) only through an overall pre-factor e^'', under the constant shift C — > C + c, (3.12) 
transform as 

C^e^^C. (3.16) 

Based on the same analysis in the previous section, the system (3.12) has an exact solution C, = const. This 
proves the conservation of C, on large scales in the generaUzed Galileon model (2.54). 



4 Conclusion 

The existence of a conserved curvature perturbation on super-Hubble scales is of particular importance in 
cosmology, which allows us relate the cun^ent observations with the theoretical predictions in the primordial 
era. The previous understanding of the conservation highly relies on the analysis of A;-essence or perfect fluid 
models. Concerning the recently proposed "generalized Galileon" models [19], it is interesting to verify if this 
conservation law still holds. 

In this work, we prove the existence of a fully nonlinear conserved curvature perturbation on large scales 
in Galileon scalar field models in two approaches. In the first approach, we show that the fluid corresponding 
to the Galileon field becomes perfect and barotropic on large scales, which are responsible to the conservation. 
The difference from fe-essence model is that, besides the energy-momentum conservation, the Einstein equation 
must be employed in order to complete the proof of barotropy. As a more familiar and simpler formalism, we 
derive the fully non-perturbative action for the curvature perturbation C, in Galileon model on large scales, and 
show that C, = const is indeed an exact solution on lai^ge scales. 

No need to say, both the two different approaches to prove the conservation are physically equivalent. The 
physical explanation of the conservation is transparent: the Galileon models differ from A;-essence model only 
with higher-order derivative terms. The effects of these higher-order derivative terms, however, vanish on large 
scales. This fact implies that, the Galileon field behaves like a A;-essencial field on large scales. Thus, on large 
scales, the Galileon field must share the same properties with fc-essencial field, including the consei^vation of 
curvature perturbation. The aim of this work, is to show this conservation in a "crystal clear" manner. 
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A A general discussion on the entropy perturbation and barotropy 

In this appendix, we will make a slightly general discussion on the barotropy of the fluid corresponding to 
single scalar field model (see [27, 32] for related discussion in A;-essence models). Especially, we will show 
that if the energy density and the pressure can be written as a system of "2-2" mapping: 

P = p{<P,4>,A(^p)), P = p{4>,4>,A(^p)), (A.l) 

where A(^p-^ and A^p^ formally represent terms which are irrelevant for our analysis in the specific problems, for 
our purpose there are higher order in spatial derivatives and can be neglected on large scales, the corresponding 
fluid always becomes barotrpic on large scales. Equivalently, the scalar perturbation becomes adiabatic on large 
scales. 

Without loss of generality, in the following we neglecte higher-order spatial derivative terms and consider 
p = p{(p, (j)) and p = p{(j), (p). Perturbing (A.l) yields 

DaP = P^Da4> + P^DaCp, DaP = P^Da4> + P4,Da4>- (A.2) 

By solving Da4> in teniis of DaP and Da4>, we get 



DaP = ciDaP + P^{^-^] Da^, (A3) 
\P'P P^J 

where we define 

= ^, (A.4) 

which has the meaning as "propagating speed of scalar perturbation". 
We can make a further variable transformation by introducing 

ea = DaP - ^Da^, (A.5) 

(A.5) can be explained as "comoving density perturbation" in single scalar field model, since which becomes 
DaP in comoving frame of scalar field where Da4> = 0. can also be explained as the "relative entropy 
perturbation" between p and (j). At this point, it is interesting to note if we alternatively introduce a "relative 
entropy perturbation" between (p and (p defined as 

Sa{4>,<P) = ^- ^, (A.6) 

(j) (p 

it is easy to show that 

€a = P^CpSa- (A.7) 

(A.5) can be used to replace DaCp in (A.3) and get 

DaP = clDaP + {cl - cl) ea, (A.8) 
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where we define 

cl = i- (A.9) 

which has the meaning as "adiabatic speed of sound". In terms of the non-adiabatic pressure, (A.8) can be 
written as an equivalent and more convenient form: 

r, = {cl - cl) ea. (A.IO) 

Thus we get the very useful conclusion: if both the energy density and the pressure can be written as functions 
of (j) and (/> as (A.l), the following three quantities are proportional to each other: 

TaCCeaCcSa- (A.ll) 

Now we are at the position to see under which condition the "fluid" becomes barotropic, i.e. DaP oc DaP 
or equivalently Ta = 0. Typically, according to (A.IO), we are left with two possibilities in order to make T^ 
vanish: 



One possibility is to make 
or equivalently 



cl = cl, (A.12) 



i = PjL = PjL. (A.13) 

P Ps 



If (A.12) or (A.13) is satisfied, the fluid becomes barotropic on all scales. In general (A.12) or (A.13) 
put constraints on the functional forms of p = p{(t), (p) and p = p{(p, (/>), and thus the structure of the 
scalar field theory. See [27, 28, 32] for a recent discussion on the barotropy of the fluid corresponding to 
A;-essence model and the constraint of exact barotropy on the functional form of P{X, (p). 

• The other possibiUty is to make 

ea = 0. (A. 14) 

As is well-known in fc-essence model, this indeed happens on large scales, where the generalized Poisson 
equation yiedls oc L'a(V^<I>) ^ 0, which yields the well-known conclusion that scalar perturbation of 
perfect fluid or fc-essence model becomes adiabatic on large scales. As we will show in this work, this is 
still the case for Galileon model, although not as explicit as /c-essence model [31-33]. See Appendix B 
for a more general discussion on the suppression of entropy perturbation on large scales. Note on large 
scales ea = implies (p = which can be viewed as the generalization of "slow-roll" condition. 

In fact, we can consider the system which can be viewed as an "N-n" mapping: 

Ql = Qliqi), i = l,---n, I=l,---N. (A.15) 

The relative entropy perturbation between any two Qi's is given by 

c fr, r, \ \^ (dQidQj dQjdQi\ 

which are linear combination of n(n — 1) /2 relative entropy perturbation among g^'s. For our purpose, one may 
find the following "N-2" mapping result be useful in practical calculations when {q,i} = {(j), (p}, the entropy 
perturbations between any of two Q/'s are 

SaiQ.Qj) = - Sa U,A . (A.17) 
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In this case, it is not surprise that all the relative entropy perturbations Sa{Qi, Qj) are proportional to each 
other This is just what happens in our system. The key point is that, in single scalar field model, all other 
entropy modes are relative entropy modes among terms with higher-order spatial gradients and thus can be 
neglected on large scales by definition, except the relative entropy perturbation between (j) and (f): Sa{4>,(t>) 
(although which itself is also suppressed by dynamical reason). Thus, Fq, ea and or any other quantity which 
has the meaning of relative entropy perturbation are just different but equivalent characterization of the same 
entropy perturbation of the system, and must be proportional to each other 

B Suppression of entropy perturbation on large scales 

In this section, we will show that, as long as the energy density and the pressure of the scalar field model can be 
written in the form p = p{(f), (j)) and P = P{4>, (j)), the entropy/non-adiabatic perturbation in the con^esponding 
system is suppressed on large scales. 

To this end, we make use of the energy and momentum constraint: 

n'^v'Gab ^\®^ + \ (^^'^ - ^abT"') = P, (B.l) 

v'Gi^hl = D'aba - lOaO = -qa, (B.2) 

combination of which yields 

^Da - ^Tafc^"') + QD'aba = DaP - Qqa- (B.3) 

(B.3) is the covariant generalization of Poission equation. The left-hand-side of (B.3) is at least third-order 
in spatial derivatives and thus is suppressed on lai^ge scales in an almost FRW background, which implies the 
right-hand-side 

la = DaP - eqa, (B.4) 

is suppressed on large scales. The next task is to show is indeed an entropy perturbation, i.e. the charac- 
terization of the non-adiabatic perturbation in our system. To see this, we have to make use of the fact that, 
the energy-flow of the scalar field theory (fc-essence, Galileons) arises only through spatial derivatives, i.e. its 
general structure is 

qa = ^Da(t>+^4^-^ + ^(^^)^ 

where q^ and q^ are general coefficients depending on specific theories, whose explicit expressions are irrele- 
vant to our following analysis. Another important observation is, qa and p, P must have the following constraint 
(neglecting higher order spatial derivative temis) 

-q^-q^- p + P, (B.6) 

which, although not obvious, can be verified explicitly through (2.61), (2.62), (2.63) and (E.l), or simply taking 
into account the functional form of the energy-momentum tensor^^ . 
Having the above relations in hand, we have 

ia -ea = - e {'^DaCl^ + "^Da^ + 0{D^)) 

4> \ (p (p J 

^eq^Sa[cl)A)+OiD^), 

Essentially, this is because the energy-momentum tensor of the scalar field is the variational derivative of a scalar Lagrangian with 
respect to the metric. 
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where we used the energy-momentum constraint p + Q{p + P) = V where V is defined in (2.9). The comoving 
density perturbation = pSa (p, (j)) is just the entropy perturbation between p and 0, which is proportional to 
Sa{4>, 0) on large scales since p = p{<j), (j)). Thus, we can conclude 

iaCC Sa{4>,(t>): (B.7) 

i.e. is the characterization of the entropy perturbation in our system. As we emphasized before, the fact that 
our system is controlled (on large scales) by two degrees of freedom and cj) implies, all different entropy 
perturbations are proportional to each other, and are suppressed. 



C Explicit expressions for ©(p), P(p) and V 



i4>) 



The discussion of this work relies on the large-scale approximation. For completeness, here we collect various 
higher-order spatial derivative teniis. 



(p) 



G 



X 



(C.l) 



and 



V 



-G 



X 



+ [-B-G^4,] Da<i)D^^ 



26 



(C.2) 



(C.3) 



D Relations among Va,ea, ia 

As have discussed in Appendix A, there are three quantities which are of particular importance in our analysis: 
the non-adiabatic pressure r^, and two covariant generalizations of comoving density perturbation ea and ia- 
In a perfect fluid or /c-essence system, they are related by 



c„ e„ 



and 



(D.l) 



In this appendix, we will derive the corresponding relations among the above three quantities in Galileon 
model (2.15). We will frequently use the definition of relative entropy perturbation between any two quantities 
SaiX,Y) defined in (2.41). 

To this end, it is useful to express them in terms of a basic quantity Sa{4>, (/)), which is the basic character- 
ization of the relative entropy perturbation in our system. We have derived the expression for p and P in terms 
of (p and (f) in (2.35)-(2.40), after some manipulations, we get: 



1 



3p,e 
26 



P L 



C<j)4>Sai<j),^)+Sa 



with 



C = P,4>P.6 - P.d>P,<t> + P.i {pM.d> - P. 6^ 



(D.l) 



(D.3) 
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and 



A 



(4>V 



-^ip)PAp) 
3p,e 



P,4>M^A4>)^'^ (%)'^(p)) +papAp)^- (An^^P)) 



26 



■7^ 



(n,A^;^^J+PA,,,A^P)Sa (7^,^(P))J • (D.4) 

5a is the summation of vaiious relative entropy perturbations among terms which are higher-order in spatial 
derivatives, and can be completely neglected on large scales. The comoving density perturbation defined in 
(2.43) is given by 



1 



3p,e 
26 



-1 



P^^(l>Sa [4>,4>) 



with 



V. 



3p,e 
26 



-1 



3p,e 
26 



nSa (Tz, 



ia defined in (2.44) is given by 



3p,e 
26 



Pi + eG,x^ 



(D.5) 
(D.6) 

(D.7) 



with 



+ 



2 

P - I ( BDa(l)D''4> + 2G,xDa<pD''X 



(D.8) 



Again, both and are higher order in spatial gradients, and should be neglected on large scales. 
E Explicit expression for t'^ 

The energy-momentum tensor corresponding to C^'^'> takes the following form: 

+CioV,Vd(/.V'^V,(/.V^Vfe(A + rff , 
where Ci, • • • , Cio are scalar coefficients given by: 

Ci = 6qJ^^ [- {acpf + {VaVbckf + 2GabV>vV - Rx] - egJJi (v.x v> - 2xa^) 

-6i^2^ fx f (VaVb0) 2 - (ncpf) + 2VaXV"X + 2□0VaXV^ 



(E.l) 



2n<^i?a6V"(/)V'' 



(3) 
X 



, (E.2) 
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Co = 6K 



(3) 



(3) 



(□</.)^ - 3 (VaVb</.)2 00 + 2 (VaVfc^)^ 



C4 
Cio 



{Ucj^Y - (VaVf,<^)" + 24^2^ 00 + 12Q 



.(3) 



12Q(Jn0-12Q^^iX-6K(3) 



(3) 



(VaVfc 

(3) 



(□0)2 + i^afcV^V* 



12K 



(3) 
X ' 



12i<:(3)n0 + 12iC^]^X - 12Q^p - ei^^^VaXV^, 



(3) 



'(3), 



-12K 



(3) 



12K^^\ 



(3) 

and represents terms whose a, fe-indices explicitly depend on or coupled to Riemannian tensors: 



^(3) 
'ab 



In deriving the above expressions, we frequently used 
as well as the Bianchi identity V aRdebc + ^bRdeca + ^ cRdeab = 0. 



(E.3) 
(E.4) 
(E.5) 

(E.6) 
(E.7) 
(E.8) 
(E.9) 
(E.IO) 
(E.ll) 



(E.12) 



F 3+1 decomposition of Riemannian quantities 

As we have seen in the bulk of this paper, besides the energy-momentum conservation equation, the Einstein 
equation must also be employed in order to complete the proof of conservation. Moreover, we frequently 
need to re-express the gravitational quantities in the equations in terms of fluid quantities. To this end, the 
"3-1-1 decompositions" of the Einstein/Riemann/Ricci tensors, which are the counterparts of the coiTesponding 
decomposition of energy-momentum tensor Tab (2.1) and VfeUa (2.7), are needed. 

These decompositions are the so-called Gauss/Codazzi/Ricci equations etc. One should keep in mind that 
they are purely kinetic equations, which are in^elevant to the underlying theory of gravitation. Some relevant 
equations which are used in this paper are summarized in the following'"*: 



Contracted Gauss eqution 



KhiR.d + hihWW'Ra'cb', 



2 1 

R-ab + g© hab + -QcFab - <7ac<^b^ ) 



(El) 



where "^Rab is the spatial Ricci tensor constructed using hab- 



'■*Here we have expressed them in terms of S, S, a°, hab, o'ab etc for convenience. For their general and more compact expressions, 
see (e.g.) [29]. 
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• Contracted Codazzi equation 

h'XRbc = Dbul - ^D.e. (F.2) 



• Ricci equation 

h'^h'^U%'^Ra'cb'd = (©^ + 39^ hab + aaUb + DaUb + CTacCrt - -^£vcrab, (F-3) 

where the last term is the Lie derivative with respect to v°- = Nu"-, with the lapse scalar. Under the 
large-scale approximation, only the first term on the right- hand-side in (F.3) survives. 

• Decomposition of Ricci scalar 

R = ^R + -Q'^ + 2Q + aabfj"^ - ^aau" - 2D„a'^, (F.4) 
3 

where ^i? is the spatial Ricci scalar. 
Some other useful relations can be derived by combining (F.1)-(F.4), e.g.: 

u'^u'^Rab = -^e^-e- aab^r"'' + aaa" + D^a'', (F.5) 

h^hfRcd = ^Rab + ^ (©^ + ©) Kb - aaab - DaQb + ^Qaab - ^aacCTb + ^£vCrab- (F.6) 

Note in the above expressions coab = is assumed, which is the case when is hypersurface orthogonal. 
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